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ABSTRACT 

We present the results from the X-ray spectral analysis of high redshift active galactic 
nuclei (AGN) in the Chandra Deep Field-South (CDFS), making use of the new 4Ms 
data set and new X-ray spectral models from Brightman & Nandra, which account for 
Compton scattering and the geometry of the circumnuclear material. Our goals are 
to ascertain to what extent the torus paradigm of local AGN is applicable at earlier 
epochs and to evaluate the evolution of the Compton thick fraction with redshift, 
important for X-ray background synthesis models and understanding the accretion 
history of the universe. In addition to the torus models, we measure the fraction 
of scattered nuclear light, / sca ttj known to be dependant on covering factor of the 
circumnuclear materal, and use this to aid in our understanding of its geometry. We 
find that the covering factor of the circumnuclear material is correlated with the line 
of sight Ah, and as such the most heavily obscured AGN are in fact also the most 
geometrically buried. We come to these conclusions from the result that / sca tt decreases 
as Ajj increases and from the prevalence of the torus model with the smallest opening 
angle as best fit model in the fits to the most obscured AGN. We find that a significant 
fraction of sources (~ 20%) in the CDFS are likely to be buried in material with close 
to 4-7T coverage having been best fit by the torus model with a 0° opening angle. 
Furthermore, we find 41 Compton thick sources in the CDFS using the new torus 
models, 29 of which we report here for the first time. We bin our sample by redshift 
in order to investigate the evolution of the Compton thick fraction by epoch. We take 
into account the incompleteness and contamination rates in the spectral identification 
of Compton thick AGN using data from simulations. We also account for the bias 
against the X-ray selection of heavily obscured sources due to flux suppression in the 
Chandra band, by restricting to intrinsic luminosities to which the CDFS is sensitive 
to Compton thick AGN (Lx~ 10 44 at z=2.5). We find a significant increase in the 
intrinsic Compton thick fraction, normalised to Lx= 10 435 ergs -1 , from « 20% in 
the local universe to «40% at z=l-4. 
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1 INTRODUCTION 

It is well known that a large fraction of AGN are obscured 
by gas and dust. This has been seen both observationally by 
the measurement of line of sight absorption in X-ray spec- 
tra (Awaki et al. 1991 Risaliti et al. 1999), and has also 



been inferred, due to a large population of obscured AGN 



being required to fit the cosmic X-ray background ( Comastri 
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et~aLlp95l |Ueda et al.||2003] |Gilli et al.||2007| ). The obscu- 
ration is, for the most part, attributable to circumnuclear 
material surrounding the black hole on parsec scales. In the 
local universe, a unification paradigm exists, which is in- 
voked to explain the differences between type 1 AGN, in 
which broad emission lines, emitted from gas orbiting close 
to the black hole, can be seen in the optical spectra, and 
type 2 AGN, in which the broad lines are not directly visible, 
but are often visible in scattered emission ( Antonucci] 1993 



Urry fc Padovani|1995 Tran|2003[). This unification scheme 
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describes a torus-like structure, which obscures the central 
engine for some sight lines, and not for others, producing 
the two observed AGN types, depending on the orientation 
of the observer. The torus paradigm of the unified scheme 
has much observational support in the local universe, but 
the relevance of the unified scheme at higher redshifts is 
not so clear ( Draper fc BallantynepOll Page et al.||2011 1. 
Characterising AGN obscuration, especially at intermediate 
redshifts and the most extremely obscured, Compton thick 
(Afa^ 1.5 x 10 24 cm~ 2 ) AGN population is key in under- 
standing the accretion history of the universe (Fabian & 
Iwasawa|1999 1 



Much can be learned about AGN obscuration through 
X-ray spectroscopy. Not only can X-rays penetrate all but 
the most severe obscuration, but the signature of obscura- 
tion is imprinted on the X-ray spectrum via soft X-ray at- 
tenuation and the emission of the iron Kq fluorescence line. 
Work to characterise obscuration for high redshift sources 
through X-ray spectral analysis has been carried out by 



Tozzi et al. ( 2006 1 using the lMs Chandra Deep Field South 
(CDFS) data and IGeorgantopoulos et al. ( 2009 1 using the 



2Ms Chandra Deep Field North (CDFN) data. These two 
fields provide the deepest pictures of the X-ray universe, 
with the latest observations in the CDFS bringing the total 
exposure in the field to 4Ms, thereby reaching on axis sen- 
sitivities of ~ 10" 17 ergs" 1 (|Xue et al.poTlf . Both 



Tozzi 



et al. (2006) and Georgantopoulos et al. ( |2009[ ) present sets 
of Compton thick candidates from their analysis, but re- 
cently more fiduciary identification of a hand full of Comp- 
ton thick sources in the CDFS have been presented by |Co-| 



mastri et al. (20111 using deep XMM-Newton data, as well 



as by Feruglio et al. (20111 and Gilli et al. (2011) using the 



4Ms Chandra data. 

In X-ray spectral analysis it is commonplace to use 
models which attribute X-ray absorption to photoelectric 
absorption along the line of sight (e.g. zwabs in xspec, see 



Mainieri et al. (2007) for an example of this.). This approach 



is sufficient at low column densities, however, when the line 
of sight column density (Ah) reaches ~ 10 23 cm -2 , the op- 
tical depth to Compton scattering becomes non-negligible. 
Not only must attenuation of the X-rays by Compton scat- 
tering be accounted for at this point, but the geometry of the 
obscuring material also plays an important role in the shape 
of the observed X-ray spectrum. For example, in Compton 
thick sources, the direct transmitted emission may not be 
visible at all, but instead scattered nuclear emission can be 
seen, either by Compton reflection from the optically thick 
wall of the torus, or by Thompson scattering from hot gas 
that fills the torus. By considering these scattering effects, 
not only can X-ray spectra indicate the level of absorption 
occurring, but they can also give clues about the geometry 
of the circumnuclear material. 

Recently, several works have presented new X-ray spec- 
tral models which properly account for Compton scatter- 



ing and the geometry of the obscuring material (Murphy 



fc Yaq oob 2009] |Ikeda et al."1|2009l |Brightman fc Nandra 



2011a). These models make use of Monte-Carlo methods to 



simulate the non-linear nature of the scattering, and include 



self consistent iron Ka emission. Murphy & Yaqoob (20091 



present a model which considers a torus with a circular cross 
section, with the opening angle fixed at 60°, whereas Ikeda 
et al. ( 2009 I present a model for a spherical torus (a sphere 



with a biconical structure removed) , with a variable opening 
angle, but where the line of sight Nn varies with the view- 
ing angle. The model we use for this work is that presented 
by Brightman & Nandra (2011al, henceforth BN11, being a 



spherical torus with a varying opening angle, but a constant 
line of sight JVh • 



Some progress has been made at using X-ray spectral 
information to constrain the geometry of the obscuring ma- 
terial in AGN. For example, 



et al. (2011) and Tazaki et al. ( 



Awaki et al. (20091, Eguchi 



2011) have applied the model 



of Ikeda et al. ( 2009k to several broad band Suzaku spectra of 



absorbed AGN, managing to place constraints on the torus 



geometrical parameters. Furthermore jRivers et al.| ( 2011 1 



have used the Murphy & Yaqoob ( 2009 1 model to constrain 



the Compton thickness of the circumnuclear material from 
an unobscured AGN, also using Suzaku. It is not only efforts 
using these new spectral models that have helped determine 
the geometry of the circum-nuclear material in AGN us- 



ing X-ray spectral data. Ueda et al. (20071 announced the 



discovery of a new type of buried AGN in which the circum- 
nuclear material is likely to be geometrically thick with a 
high covering fraction. This conclusion was reached by the 
measurement of a very small (< 0.5%) scattered fraction, 
/scatt- This so called "scattered" emission is thought to orig- 
inate from the Thompson scattering of the primary X-ray 
photons by hot electrons within the cone of the torus, and is 
ubiquitous in obscured AGN ( Turner et al.||1997 l, or emis- 
sion from circum-nuclear plasma photo-ionised by the AGN 



(Guainazzi & Bianchi|2007l. A small scattered fraction im- 



plies a small opening angle of the torus, or otherwise an 
under abundance of the gas responsible for the scattering. 
This picture is supported by Noguchi et al.] | |2010[ ) who found 
that the scattered fraction correlates with the [O ill] narrow 
emission line to X-ray flux ratio. As the [O ill] emission line 
also originates in gas responsible for the scattering of the 
X-rays, it is also gives an indication of the covering factor 
of the gas, or its abundance. 



While these works have concentrated on a small set of 
high quality X-ray spectra of local AGN, the aim of this pa- 
per is to take advantage of the 4Ms CDFS data, which has 
resulted in many spectra containing of order of hundreds of 
counts or more, enabling application of more such sophisti- 
cated models. Our goal is to use both these new X-ray spec- 
tral models and the measurement of /scan to characterise the 
nature of the obscuring material in AGN in the CDFS, and 
to use them to identify Compton thick sources. In section [2] 
we describe our data set and the data reduction techniques, 
in section [3] we describe our spectral fitting technique using 
the new torus models of BN11, and in section|4]we present 
our results on our investigation into the geometry of the ab- 
sorber in high redshift AGN and identification of Compton 
thick sources. We discuss these results in section In this 
work we assume a flat cosmol ogical model with _Hp= 70.5 



km s _1 Mpc" 1 and fi A =0.73 (Komatsu et al. 20091. For 



measurement uncertainties on our spectral fit parameters 
we present the 90% confidence limits given two interesting 
parameters. 



2 CHANDRA DATA 

The Chandra Deep Field South currently consists of 52 in- 
dividual observations, giving ~4 Ms of total exposure time 
in the field. Of these observations, nine were made in 2000 



(~1 Ms,|Giacconi et al.|2002) , 12 made in 2007 (~ 1 Ms,|Luo 
et al.|2008[ ) and a further 31 made in 2010 (~2 Ms, |Xue et ah 
2011[ ). We analyse the entire 4Ms data set, where data were 
screened for hot pixels and cosmic afterglows as described 



in Laird et al. (20091. For our spectral analysis, we use the 
source list presented by |Luo et aL| ( |2008[ ) , based on detec- 
tions in the 2Ms data set, consisting of 462 sources down to 
a flux limit of ~ 1 X 10" 16 — - — -a " _1 
band and ~ 2 x 10~ 17 ergs cm 
As the focus of this work in on X-ray spectral analysis, we 
neglect the faintest sources newly detected in the 4Ms data 



ergs cm " s * in the 2-8 keV 
2 s" 1 in the 0.5-2 keV band. 



set presented by Xue et al. (20111. Events were extracted 



from t he data using the ACTS extract (AE) software pack- 
age rl (|Brooseti^l2010 1 . AE extracts spectral information 



for each source from each individual observation ID (obsID) 
based on the shape of the local point spread function (PSF) 
for that particular position on the detector. We choose to 
use regions where 90% of the PSF has been enclosed at 1.5 
keV. Background spectra are extracted from an events list 
which has been masked of all 776 detected point sources in 



Xue et al. (20111, using regions which contain at least 100 



counts. AE also constructs response matrix files (RMF) and 
auxiliary matrix files (ARF). The data from each obsID is 
then merged to create a single source spectrum, background 
spectrum, RMF and ARF for each source. 



Of the 462 sources presented by Luo et al. ( 2008 1, we use 



408 secure spectroscopic or photometric redshifts as com- 



piled by Silverman et al. ( 2010 1 using new data from the 



VLT and Keck as well as archival data from|Szokoly et"aL 
( 2004| , [LeFevre et al.|||2004||,|Wolf et al.|(|2004|>,|Zheng et al 



(20041, Ravikumar et al. (20071, Wolf et al. (20081, Popesso 



et al.| ( |2009| >, |Treister et al.| ( |2009[ ) and |Balestra et al.| ( |2010 |. 
We use a further 41 photometric redshifts presented in|Xue 



et al. 



(20111 for sources not presented in Silverman et al 



(20101. The total number of sources for which we have ob- 
tained a redshift for is 449, of which 296 are spectroscopic 
redshifts. Considering both spectroscopic and photometric 
redshifts, the redshift identification completeness of our sam- 
ple is 97%. For the remaining 13 sources without redshifts, 
we neglect them from our analysis. Fig. IT] shows the distri- 
bution of spectral counts for our sample, and Fig. [2] shows 
their redshift distribution. 



3 X-RAY SPECTRAL ANALYSIS 

3.1 X-ray spectral fitting 

We carefully assess our options for our spectral fitting tech- 
nique. Due to relatively low count statistics in the CDFS 
spectra, and our wish to preserve as much spectral resolution 
as possible, we do not wish to significantly group our spec- 
tra. The standard fit statistic used in X-ray spectral analysis 
is the \ 2 statistic. However, this is known to fail when the 



1 The ACIS EXTRACT software package and User's Guide are avail- 
able at http://www.astro.psu.edu/xray/acis/acis_analysis.html 
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Figure 1. A histogram of the total number of 0.5-8 keV spectral 
data counts in our sample. 
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Figure 2. Histograms of the redshifts of the sources in our sam- 
ple, shown with both linear and logarithmic binning. The hatched 
histograms within the empty histograms show the photometric 
redshifts used. 



spectrum contains bins in which there are only a few counts, 
making it necessary to bin the spectra, usually with a mini- 
mum of 20 counts per bin. The fit statistic most appropriate 
for fitting with low numbers of counts per bin is the Cash 
statistic (cstat, Cash|[l979 1 which uses a Poisson likelihood 
function. For the spectral fitting in the CDFS carried out 



by Tozzi et al. (20061, the authors carried out simulations 



of low count spectra to compare results from fitting using 
the Cash statistic and the \ 2 statistic. They find that not 
only does the Cash statistic produce a result closer to the in- 
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put parameters, but the rms dispersion on the result is also 
lower when using Cash statistics. They also emphasise the 
point that as the \ 2 statistic requires that the spectrum be 
binned, the spectrum used has essentially been smoothed re- 
moving fine spectra details, such as emission lines. The Cash 
statistic is not normally suitable for background subtracted 
spectra, nevertheless, xspec has a modified version of cstat 
for when a background spectrum is included 2 ! We use xspec 
version 12.6.0q to carry out X-ray spectral fitting. For these 
reasons, we have selected to use the Cash statistic in our 
analysis, where we only lightly group our spectra with at 
least one count per bin as suggested in the documentation, 
using the HEASARC tool grppha. 



3.2 Determining the geometry of the absorber 

A major motivation behind this work is to assess the geom- 
etry of the obscuring material in AGN in the CDFS using 
X-ray spectroscopy. For this we use the new torus models 
described in BN11. These models are based on Monte-Carlo 
simulations of an X-ray point source, emitting a power-law 
spectrum, Fe — AE~ r (range 1.0 ^ T «J 3.0) at the centre of 
a distribution of a cold neutral medium, with a line of sight 
column density to the observer, _Vh (range 10 20 ^-/VhSJ 10 26 
cm -2 ). The interactions modelled are photoelectric absorp- 
tion, Compton scattering, and K-shell fluorescence from sev- 
eral elements, including Fe Ka and Fe K/3 emission, where 
the abundances of the elemental constituents are assumed 
to be solar. The geometry of the obscuring material mod- 
elled is a spherical torus which is essentially a sphere, with 
a biconical cone cut into it. The geometrical parameters 
of this distribution are the half opening angle of the torus 
(25.8° ^ 9 t0 r ^ 84.3°) and the inclination angle of the torus 
(18.2° < 9i sC 87.1°). Fig. [3] shows the geometry of the 
torus, and shows how 8t or and 9i are defined. For sight lines 
through the torus, the _Vh is independent of the viewing an- 
gle. Also presented in BN11 is a special case, whereby the 
source is covered by 4-7T steradians of material in a spherical 
distribution. 

It is unlikely however, that the torus geometrical param- 
eters can be constrained for many individual faint sources 
in the CDFS. Therefore, in our spectral analysis, our goal 
is to investigate the prevalence of the best fit model, con- 
sidering three geometrical scenarios, where the geometrical 
parameters have been fixed. These are a torus geometry, 
where the opening angle is fixed to either 60°, 30° or 0° 
(i.e spherical) and the inclination angle is fixed to edge on. 
For determining the geometry of the absorber, we fit only 
the 408 sources with the most secure redshifts, as presented 
in [Silverman et al.| ( |2010[ ) so as to minimise the effect from 
uncertainties in the photometric redshifts. We fit each of 
these source spectra with the following five combinations of 
spectral models, which are physically motivated and rep- 
resent different geometrical scenarios for the circumnuclear 
material. For the torus scenarios where the opening angle 
is >0°, we include a secondary power-law in the fit, which 
represents a fraction of the nuclear emission posited to be 
scattered into the line of sight by gas filling the cone of the 
torus. 



https://astrophysics.gsfc.nasa.gov/XSPECwiki/low_count_spectra 





Figure 3. Geometry of the torus model presented in BN11. 
The shading only serves to differentiate the different opening an- 
gles used, and does not correspond to density. Typical spectra 
are shown below the diagram. The first panel shows the unob- 
scured sight-line (6>i nc =20°), whereas the other panels show ob- 
scured sight-lines (_Vh=10 23 
B inc =m° and Ar H =10 25 cm- 
black lines are for 8tor=30 c 

6»tor=60° 



cm~ 2 , inc =8O°; 7V H =10 24 cm" 2 , 

2 , 0inc=8O° respectively). The solid 

and the dashed red lines are for 



• A) The torus model of BN11 with a fixed opening an- 
gle of 60° , and edge on orientation, accompanied by a scat- 
tered power-law, with r scat t fixed to the value of the primary 
power-law. This model has four free parameters, _Vh; F; and 
the normalisations of the power-laws, Ai and A2 

• B) The torus model of BN11 with a fixed opening angle 
of 30° , and edge on orientation, accompanied by a scattered 
power-law, with V fixed to the value of the primary power- 
law. This model has four free parameters, _Vh; T; Ai and 
A 2 

• C) The spherical model of BN11, a scenario in which 
the X-ray source is completely covered with 47r steradians of 
obscuring material. This model has three free parameters, 
_Vh; r and Ai. We include no scattered component for this 
model as it represents the case where there is no escape 
route for the primary radiation to be scattered into the line 
of sight. We do however calculate an upper limit on this 
component. 

The new torus models take the place of well established 
models used previously to model absorption (e.g. zwabs) and 



reflection from Compton thick material (e.g. pexrav). The 
model plcabs ( Yaqoob|1997 l describes the special case where 
the source is covered by 4n steradians of obscuring material, 
analogous to our spherical model. 

The final two models describe unobscured X-ray emission, 
though these sources are not the focus of our work. 

• D) simple power-law model, with two free parameters, 
the power-law index, F and the normalisation, Ai, repre- 
senting unobscured X-ray emission. 

• E) The torus model of BN11 with a fixed opening angle 
of 60°, and face on orientation. This model is essentially a 
simple power-law, but includes the reflection features from 
the torus, such as the Fe Ka line and Compton hump. This 
model has three free parameters, Nh (through the torus); V 
and Ai 

We choose the model which presents the best c-stat 
value, after accounting for the number of free parameters 
of the model. For example, for model C, with three free pa- 
rameters to be chosen as best fit model over model D, with 
two free parameters, it must present a c-stat improvement 
of at least 1.0, where 1.0 is the Ac-stat associated with the 
addition of a free parameter at the 1-er confidence level, i.e. 
cstat(C)<cstat(D)-1.0 

For spectra with a low number of counts, it can be dif- 
ficult to obtain a converged spectral fit due to large data 
uncertainties. We perform error analyses on each model fit, 
using the xspec error command, for two interesting pa- 
rameters and Ac-stat=4.61, which corresponds to 90% con- 
fidence range. Despite these actions, it is still difficult to 
constrain the fit parameters well with low count statistics. 
Fig H] shows the binned average uncertainty on V as a func- 
tion of total 0.5-8 keV source counts. Below ~ 600 counts, 
< Ar >> 0.6, which could lead to higher uncertainties in 
the other parameters. For this reason, we choose to fix V 



at 1.9, the canonical value for Seyfert galaxies (Nandra & 
Pounds||l994 |, for spectra with less than 600 counts. How- 
ever, in this case, if the best fit model for which F remains 
free provides a significantly better fit at the 90% confidence 
level (Ac-stat> 2.71) than for the best fit model where T 
is fixed at 1.9, we select the best fit model where F re- 
mains free. The motivation behind including this additional 
criterion for allowing T to be free when there is less than 
600 counts is to account for sources with extreme F values. 
These include for example, Blazars, which can present very 
high (~ 3) to very low (~ 1) F values ( Donato et al.| 2005). 
Without this consideration, these sources may be identified 
as Compton thick if F were fixed to 1.9, due to their flat 
spectra. This is done for 60/310 sources with less than 600 
counts. However, more than half of these have more than 
200 counts, where V can still be fairly well constrained. 



3.3 Identifying Compton thick sources 

The second major aim of this work is to identify Compton 
thick sources in the CDFS using these new spectral models. 
Strictly a Compton thick source is that with Nn> 1.5 x 10 24 
cm" 2 , however, for simplicity here we use 10 24 cm" 2 as the 
lower limit. We then aim to assess the Compton thick frac- 
tion of sources (sources with 10 24 <Nn< 10 26 cm -2 with 
respect to those with 10 20 <N H < 10 26 cm" 2 ) at high red- 
shift. Spectral fitting with models which correctly account 
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Figure 4. The average uncertainty on T as a function of the 
number of spectral data counts for the sources in our sample. The 
dotted line marks 600 counts, the minimum number of counts we 
require to leave T free for fitting 



for Compton scattering and the geometry of the absorber 
and which also include self consistent Fe Ka emission are 
valuable for identifying Compton thick sources, as the equiv- 
alent width of the Fe Ka emission is tightly linked to Nu 



( Leahy fe Creighton|1993l|Ghisellini et al.|1994||Brightman 



fe Nandra||2011a l. Furthermore, intrinsic source luminosity 
calculations are more realistic when using these new mod- 
els. Models based on a slab geometry such as pexrav un- 
derestimate the intrinsic source luminosity by several fac- 
tors when fitting Compton thick sources (|Murphy fe Yaqoob| 



2009| |Brightman fe Nandra1|2011a[ ). This is due to the fact 



that pexrav assumes a slab of optically thick material that 
subtends a solid angle of 2-7T steradians from the perspective 
of a point source, relevant for an optically thick accretion 
disk. However, the torus is likely to subtend a much smaller 
solid angle, and thus the reflection strength is weaker with 
respect to the accretion disk. Therefore the intrinsic source 
luminosity for a torus model is higher than that for pexrav 
for the same reflection strength. For identifying Compton 
thick sources, we utilise all 449 sources with available red- 
shifts, be it spectroscopic or photometric. 

We also fit our spectra with standard analytical mod- 
els to compare to the results from the torus models. For 
this purpose we use a simple power-law model, combined 
with zwabs ( Morrison fe McCammon|1983 l to model photo- 
electric absorption, zgauss (line energy fixed to rest frame 
6.4 keV, width fixed to 0.001 keV) to model any Fe Ka 



emission and pexmon (Nandra et al. 20071 to model reflec- 



tion and associated Fe K emission. We also use a secondary 
power-law component to fit the scattered emission, where 
the power-law index is fixed to that of the primary power- 
law component. Once again, for the addition of any free 
parameters, the c-statistic must improve by 1.0. In Fig [5] 
we show a comparison of the Nh as measured by the torus 
models, with the JVh measured by the analytical models for 
the most absorbed sources. This shows that the analytical 
models systematically overestimate the Nh with respect to 
the torus models as they do not account for attenuation of 
the spectrum by Compton scattering. 
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33 




geometrically thick material with close to total coverage as 
they are best fit by the 0° torus model and do not require a 
scattered component. These sources are analogous to the so 
called "new type" AGN as discovered in the local universe 



by Ueda et al. ( 2007 1 



N H (cm" 2 ) 



Figure 5. The Nh as determined by the torus model (x-axis) 
against the Nh as determined by the analytical models, for heav- 
ily absorbed sources. It can be seen that the analytical models 
overestimate the TVh with respect to the torus models as they do 
not account for attenuation of the spectrum by Compton scatter- 
ing. 



4 RESULTS 

4.1 Geometry of the absorber 

In Table IT] we present the best fit nrodel data, which shows 
the number of spectra where each model combination has 
been selected as the best fit model. For low Ah, models 
A, B and C, which model obscured sight lines, are largely 
degenerate, as only the geometry of the absorber differenti- 
ates them. The geometry has negligible effect below ~ 10 23 
cm" 2 as Compton scattering is not important. In column 
3 we show the number of each model selected where the 
Ah is measured to be > 10 23 cm -2 . If we consider the 
fraction of scattered nuclear X-rays, / sca tt, which has been 
shown to be dependant on the opening angle of the torus, 
we find that heavily absorbed sources with a scattered frac- 
tion /scatt > 1.0% are best fit by the torus model with the 
largest (60°) opening angle (model E) over the torus model 
with the smallest (30°) opening angle. Lastly, we present 
the mean / SC att for the two torus models, for spectral fits 
where there are at least 600 counts in the spectrum. Model 
B with the 30° opening angle gives a systematically smaller 
/scatt (1.4 ± 0.3%) than model A with the 60° opening an- 
gle (2.1 ± 0.7%). The average / sca tt we measure for all fits 
that require the scattered component with > 600 counts is 
1.7 ± 0.4%. These results present us with a self consistent 
picture, whereby the fraction of scattered X-rays measured 
in the spectra is closely linked to the covering factor of the 
torus. Our results show that while the torus paradigm is still 
applicable at high redshifts, there is a significant fraction of 
sources (~ 20%) where the source is likely to be buried in 



Having found that our results are consistent with / sca tt 
being linked to the opening angle of the torus, we can in- 
vestigate how the covering factor varies with other source 
parameters as well as how the fraction of fits which prefer a 
torus with high covering factor over a lower covering factor 
behaves. In Fig.lSJ we plot / sca tt against Ah- We plot as thick 
red bars only sources with 10 43 <Lx< 10 44 ergs -1 and 
1 < z < 3 in order to break any degeneracy with these pa- 
rameters. We do not explore the relationship between / sca tt 
and Lx as there may be contamination of scattered X-rays 
from star forming processes in the host galaxy, the relative 
strength of which are likely to be linked to the luminosity 
of the AGN. Furthermore the constraint on / SC att decreases 
towards high redshift due to the scattered component being 
redshifted out of the Chandra band, and for this reason we 
do not explore the relationship between / SC att and redshift. 

Spectra where /scatt > 10% are likely to represent some 
scenario where the source is only partially covered, or where 
the Ah has changed during the 10 years of observations (e.g. 
NGC 1365, |Risaliti et "aL"1|2009| ), and thus not relevant to 
the scattered X-rays we arc considering here. We therefore 
disregard sources with / SC att> 10% in our analysis of the 
scattered fraction. We plot the binned average of /scatt for 
the various sub-samples in the lower panels. We find that 
/scatt is strongly dependent on A H , with / SC att decreasing 
with increasing Ah, implying that the covering factor of the 
torus increases with line of sight obscuration. 

To investigate the effect of spectral variability in more 
detail, we look for photometric variablity in the most ob- 
scured sources (Ah> 10 23 cm" 2 ), which is likely to be 
caused by spectral variations rather than changes in the 
intrinsic source power. During spectral extraction, the soft- 
ware package ACIS extract computes the significance of the 
Kolmogorov-Smirnov (KS) statistic when comparing a uni- 
form count rate model to the observed distribution of source 
event time stamps. For heavily absorbed sources where this 
KS statistic is < 0.1, suggestive of variability, we split the ob- 
servation into two, analysing the observations made in 2000 
and 2007 (2Ms) and the observations made in 2010 (2Ms). 
For 22 sources where the Ah has changed by AAh>1 dex, 
from the first 2Ms to the second, we remove them from our 
analysis of the scattered fraction. We find that removing 
these sources has no effect on our conclusions regarding the 
dependance of / sca tt on A H . 

We also investigate the relative prevalence of the two 
torus models with different opening angles, and how this 
varies with Ah, Lx and redshift. In Fig [7] we plot the frac- 
tion of sources for which the torus model with a 30° opening 
angle is preferred as best fit over the torus model with a 60° 
opening angle for sources with Ah> 10 23 ' 5 cm -2 , as at low 
Ah values these models are degenerate as Compton scatter- 
ing is not significant. Again, in order to break degeneracies, 
for Lx we plot only for 1 < z < 3, and against redshift we 
plot only for 10 42 <Lx< 10 44 ergs" 1 . These results also 
indicate that the most heavily obscured sources are in fact 
the most heavily buried, as the torus model with the small- 
est opening angle is preferred for the most heavily obscured 



Table 1. Details of the best fit models. Column (1) gives the model description; Column (2) gives the number where each model 
combination is the best fit; Column (3) gives the number where each model combination is the best fit and N^> 10 23 cm~ 2 ; Column (4) 
gives the number where each model combination is the best fit, Nn> 10 23 cm — 2 and /scatt is > 1%; Column (5) is the average / SC att for 
each, where the total number of counts in the spectrum is > 600 



Model 



(1) 



All 7V H > 10 23 cm" 2 Nn> 10 23 c 



(2) 



(3) 



+/s, 



(A) 60° torus 

(B) 30° torus 

(C) 0° torus (sphere) 

(D) power-law 

(E) unobscured torus 



77 (19%) 
83 (20%) 

135 (33%) 
90 (22%) 

23 (5.6%) 



61 (42%) 
56 (39%) 
28 (19%) 



Total 



408 



1 15 



>1% 
(4) 



< /scatt > (%) 

(5) 



50 (34%) 
40 (28%) 



90 (62%) 



2.1 ± 0.7% 
1.4 ±0.3% 



1.7 ±0.1 



sources, whereas the torus model with the largest opening 
angle is preferred in the least obscured sources. This sup- 
ports our conclusion from the / sca tt data. 



4.2 Identifying Compton thick sources 

Of the 449 sources we study here, we find that 41 of them 
are Compton thick from our analysis, giving an observed 
Compton thick fraction for all sources of 9%. We present 
the details of the spectral fit parameters, including error 
analyses, of 20 of these sources which have been found at 
the 90% confidence limit (Ac-stat criterion of 2.71), and are 
thus the most secure. These are presented in Table[2j and the 
spectra are presented in Figs [8] [9] and [To] These Compton 
thick sources are predominantly best fit by the torus model 
with the 30° opening angle, and have very small / SC att values 
(11 fits with /scatt < 1.0%), which supports our conclusion 
that the most heavily absorbed sources are in fact the most 
geometrically buried as well. Due to the low count nature 
of these sources, the power-law index has been frozen to 1.9 
in all but four fits. A total of nine of the Compton thick 
sources have intrinsic Lx> 10 44 ergs -1 , making them type 
2 quasars. 

The highest redshift Compton thick source we find is 
at z = 7.62 (306). The redshift is photometric and given 
by Luo et al. (20101. However, an alternative photometric 

If we allow 



redshift of 3.3 is also presented by the authors, 
the redshift in the spectral fit to be a free parameter, we find 
that the best fit redshift is in fact 6.9l 1 ' 4 , which is consistent 
with the preferred photometric redshift in |Luo et al.| | |2010l |. 
If the source were indeed at z=7.62, this would make this 
source the highest redshift Compton thick AGN to date and 
with an intrinsic deabsorbed 2-8 keV luminosity of 3.0 x 10 44 
ergs -1 . We present the spectral fit details and spectra for 
both photometric redshifts given in Luo et al. (20101. 



4.3 Evolution of the Compton thick fraction 

Previous works have shown that the fraction of absorbed 
AGN increases with redshift (e.g. Hasingcr 2008). The wide 
redshift range of the sources in the CDFS allows for an in- 
vestigation into if and how the Compton thick fraction of the 
AGN population evolves with epoch. We decide to investi- 
gate this evolution with three redshift bins, z=0.1-l,l-2.5 






OB 

o 
to 




0.02 
0.00 



10 



22 



10 



23 



10 



24 



N 



H 



Figure 6. Upper panel: The vertical bars represent the 90% con- 
fidence limit on the range of /scatt for each source, against JVh- 
Red data points are for sources with 10 43 <Lx< 10 44 ergs -1 
and 1 < z < 3 in order to break any degeneracy with those pa- 
rameters. Lower panels: The binned average of /scatt against Afjj 
for the subsample defined above (red bars in the upper panel). 
The error bars represent the standard error (i.e. the standard de- 
viation of the data divided by the square root of the number of 
data points) on the / SC att values. 
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Figure 7. The fraction of sources where the torus model with 30° opening angle is chosen as best fit model over the torus model with 
60° opening, against TVh and Lx and redshift. We have plotted only for sources with JVh> 10 23 ' 5 cm -2 , as at low Nn values these 
models are degenerate as Compton scattering is not significant. Against Lx we plot only for 1 < z < 3, and against redshift we plot only 
for 10 42 <L X < 10 44 ergs" 1 

Table 2. Details of the secure Compton thick sources found in this analysis. Column (1) is the Luo, et al (2008) source number; Column 
(2) is the number of 0.5-8 keV spectral counts; Column (3) is the redshift of the source; Column (4) is the best fit model for the spectrum 
as described in sectional Column (5) is the Cash statistic of the fit; Column (6) is the measured Nn in units of 10 24 cm -2 ; Column (7) 
is the photon index, T, where ' indicates that the index was frozen in the fit; Column (8) gives the scattered fraction, /scatt; Column (9) 
gives the unabsorbed rest-frame 2-10 keV luminosity of the source; Column (10) gives the redshift reference from Table 4 of Silverman, et 
al (2010). The redshift for source 306 was taken from Luo, et al (2010) and is a photometric redshift. For this source we have presented 
spectral fits for the photometric redshift presented (7.62), and the alternative photometric redshift given by the authors (3.3). The 
following sources are in addition to those presented below and have been detected with 1-<t significance - 8, 82, 121, 162, 170, 201, 208, 
223, 258, 260, 262, 265, 313, 314, 316, 330, 337, 339, 376, 390 and 399. 
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Figure 8. Unfolded spectra of the secure Compton thick sources found in this analysis, using the torus models of BN11. The red dashed 
lines show the torus component itself, which includes self consistent Fe Ka emission, whereas the black dashed lines show the scattered 
nuclear component, which typically accounts for < 1% of the direct emission in these sources. 
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Figure 9. Unfolded spectra of the secure Compton thick sources found in this analysis, using the torus models of BN11. The red dashed 
lines show the torus component itself, which includes self consistent Fe Ka emission, whereas the black dashed lines show the scattered 
nuclear component, which typically accounts for < 1% of the direct emission in these sources. 
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and 2.5-4, designed to cover the redshift range of the CDFS. 
In order to calculate the intrinsic Compton thick fraction 
of each bin, we must take into account several biases and 
limiting factors in our data. Firstly we must determine the 
reliability of our spectral fitting method to correctly iden- 
tify Compton thick sources, especially at the faintest fluxes, 
which we do via simulations. Secondly, as the sample is X- 
ray selected, it is biased against the most absorbed sources, 
due to severe flux suppression in the X-ray band. Further- 
more, the survey solid angle is a sharply de creasing function 
of flux below Fo.5-8 ~ 10" 15 erg cm" 2 s" 1 JLuo et al.|2008k, 



which means that it is further biased against faint absorbed 
sources, for which a correction must be made. The source 



catalogue of Luo et al. (2008) is a band merged sample, de- 



tected in the soft (0.5-2 keV), hard (2-8 keV) and full (0.5- 
8 keV) bands. To account for the factors described above, 
we require a single band selected sample. We thus use the 
full band selected sample, as it has the greatest number of 
sources (406). 

4-3.1 Simulations 

Above we have described the details of the most secure 
Compton thick sources from our analysis, at the 90% con- 
fidence limit. In our spectral analysis, we adopted the 1- 
a significance level when adding parameters to the spec- 
tral model, such as absorption. We run simulations in or- 
der to ascertain the fraction of Compton thick AGN which 
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Figure 10. Unfolded spectrum of source 306, found to be Compton thick in our analysis, given a photometric redshift of 7.62 by Luo, 
et al (2010). We have not binned the spectrum for plotting here so as to show the line at ~ 0.7 keV, which is consistent with Fe Ka 
emission at the photometric redshift. The spectrum of this source at the alternative photometric redshift given by Luo, et al (2010) of 
3.3 is also shown for comparison. 



are correctly identified using our spectral fitting method, as 
well as to determine the number of contaminating sources, 
sources which are incorrectly identified as Compton thick. 
We use the fakeit command in xspec to produce sim- 
ulated spectra, assuming a torus model with 30° open- 
ing angle and a 1% scattered component, using the Chan- 
dra response and a typical background file from a typical 
source in the CDFS. We simulate 100 spectra for each Nh 
(10 20 ' 5 , 10 215 , TO 22 ' 5 , 10 23 ' 5 , IO 245 , TO 25 ' 5 ), for three different 
0.5-8 keV count rates (10~ 4 ' 5 , 10~ 3 ' 5 , IO" 25 s _1 ) and three 
redshifts (z=0.3, 1.0 and 3.0). The rate of Compton thick 
identification and the contamination rate determined from 
these simulations are plotted in Fig. |11| as well as the distri- 
bution of count rate in three redshift ranges, corresponding 
to the simulations. We find that, as expected, the rate of 
Compton thick identification is at its lowest at < 40% at the 
faintest count rates, but also in the lowest redshift sources. 
This of course improves with increasing count rate and also 
as the redshift increases. A significant increase in the rate of 
identification is seen from redshift 0.3 to 1, due to the soft 
excess emission being redshifted out of the band pass and 
the reflected light being redshifted into it. We thus also plot 
the interpolation of the identification rate for z=0.7. The 
contamination rate is < 10% for all count rates and red- 
shifts. Plotted in Fig. [12] is the rate of CT identification for 
a source with an intrinsic _Lx=10 43 ergs -1 , as a function of 
redshift. This shows that despite the Chandra band observ- 
ing a higher rest-frame energy for higher redshift sources, 
and thus being less affected by absorption, identification of 
Compton thick AGN becomes more difficult as redshift in- 
creases, for a constant luminosity. 

While simulated data can give us statistical information 
on the rate of Compton thick identification, it is difficult 
to determine which individual sources have been correctly 
identified, or not identified at all. One test that can be per- 
formed is to check if optical broad lines have been detected 
in any of the sources identified as Compton thick, as these 
sources are unlikely to be Compton thick in reality. We find 
that 9 such sources are optical broad line sources (given in 
Silverman et al.|2010 1 and have been identified as Compton 



thick in our analysis. We therefore remove these from our list 
of Compton thick candidates. However, this contamination 
rate is consistent with data from our simulations. Given our 
knowledge about the rate of Compton thick identification 
from these simulations and the rate of contamination, we 
correct on a statistical basis, the observed number of Comp- 
ton thick sources. For the z=l-2.5 bin we use the z=l sim- 
ulation results to correct the numbers, and for the z=2.5-4 
bin, we use the z=3 results. As we have shown from our sim- 
ulations, the rate of Compton thick identification increases 
significantly from z=0.3 to z=l. The redshift distribution of 
our z=0.1-l sources peak around z=0.7, and thus to avoid 
underestimating the identification rate, we use the interpo- 
lation of the z=0.3 and z=l results to correct the numbers 
in this bin. 



4-3.2 Correcting for survey biases 

As our sample is X-ray selected, it is intrinsically biased 
against heavily obscured AGN, due to the significant flux 
suppression occurring in these sources. In order to account 
for this bias, we restrict our analysis to an intrinsic lumi- 
nosity range, where the lower limit is determined by the in- 
trinsic luminosity to which the survey is sensitive to Comp- 
ton thick sources. Below IO -5 counts s _1 in the 0.5-8 keV 
band, the survey solid angle becomes very small, which 
could introduce large uncertainties in our data due to small 
number statistics. We therefore apply a count rate cut at 
this level to our sample. Fig. [13] plots the intrinsic 2-10 
keV luminosity against redsift of the Compton thick AGN 
above this count rate limit, which has been determined from 
the spectral fit model. The solid curves on this plot shows 
the intrinsic luminosity we are sensitive to sources up to 
iV H = 10 24 ,10 24 ' 5 ,10 25 cm" 2 with this adopted count rate 
limit. We then calculate the intrinsic Compton thick frac- 
tion by restricting our analysis to sources with intrinsic lumi- 
nosities above the curve for Nh_= 10 25 cm~ 2 ,and Lx> 10 42 
ergs -1 in order not to include star forming galaxies. While 
using this curve we loose a few CTAGN in the highest red- 
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Figure 11. Top panel: Solid lines - the rate of correct identifi- 
cation of CT sources from our simulations for a range of count 
rates, for sources at z=0.3,1.0,3.0. The dot-dashed line is an in- 
terpolation of the z=0.3 and z=l results to z=0.7; dotted lines 
- the fraction of Compton thin sources incorrectly identified as 
Compton thick for the same redshifts. Bottom panels: The distri- 
bution in 0.5-8 keV count rates for sources at z=0.1-l, 1-2.5 and 
2.5-4. The filled histogram shows the Compton thick sources 
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Figure 12. The rate of correct identification of CT sources from 
our simulations for a source with Lx=10 43 ergs -1 as a function 
of redshift. 
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Figure 13. The intrinsic 2-10 keV luminosity of the CTAGN 
with a 0.5-8 keV count rate of > 10~ 5 counts s —1 in our sam- 
ple against redshift. The vertical dashed lines show the redshift 
binning we use, z=0. 1-1, 1-2. 5, 2. 5-4, while the solid curve repre- 
sents the intrinsic luminosity we are sensitive to sources up to 



10 24 10 24.5 10 25 



with the adopted count rate limit. 



shift bins, we can be certain that we are complete to all 
CTAGN. 

To correct for the differential survey solid angle as a 
function of flux, we bin up our sample into three bins of 
1 dex in 0.5-8 keV count rate, and convert the observed 
number of sources into number densities using the average 
survey solid angle for the 0.5-8 keV band within that bin. 
The solid angle data presented in |Luo et ah| ( |2008[ ) is given 
as a function of flux, which has been converted from a count 
rate with a power-law of T = 1.4, and so we convert this 
flux back into a count rate using the same model, so that 
we may carry out the integration. As mentioned above, we 
apply a count rate cut of 10 -5 counts s _1 due to the very 
small survey solid angle below this limit. 

Once we have calculated number densities, having cor- 
rected for all biases aforementioned, we can then calcu- 



late the Compton thick fraction for each redshift bin. 
For each redshift bin a different range in luminosity has 
been used, due to the increasing luminosity sensitivity to 
Compton thick AGN with redshift. However, X-ray ob- 
scuration is heavily dependent on X-ray luminosity, so in 
order to evaluate Compton thick fractions as a function 
of redshift, we must determine the Compton thick frac- 
tion at a common luminosity. This was done similarly by 



Hasinger (2008), H08 henceforth, who evaluated the de- 



pendence of AGN obscuration on luminosity and redshift. 
The author finds that fa is a linear function of logioLx 
with a slope of -0.28, where the constant increases with 
redshift. We require a relation for /ct vs. Lx which we 
convert from /2 vs. Lx • In H08, X-ray absorption was 
determined through the use of hardness ratios, such that 
/ 2 = JV(21.5 <logi Af H < 24)/AT(20 <logi Af H < 24). We as- 
sume that most of that sample is Compton thin and does 
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Figure 14. The evolution of the Compton thick fraction with 
redshift, which has been normalised at logio-Lx=43.5 in each red- 
shift bin, using data from H08. The 1-2.5 and 2.5-4 redshift bins 
have been merged here to improve the statistics. The uncertain- 
ties on /ct have been propagated from the Poisson errors on 
N(N H > 10 24 cm" 2 ) and N(N n < 10 24 cm" 2 ) Th e Swift/BAT 
sample has an average luminosity of logioLx=43.6 | |Burlon et al.| 
|2011| ) , and so we plot the data from this sample with no renor- 
malisation. 



not contain many sources with Nn> 10 24 cm 2 . In order 
to convert the relation for fo into a relation for /ct, where 
/ct = iV(24 <logi A^H< 26)/iV(20 <log 10 N H < 26), we as- 
sume that iV(24 <logi iVH< 26)=JV(21.5 <logi iVH< 24), 
which is an almost fiat Nu distribution above Nn—10 2 ' 5 
cm -2 . This gives us the relation /ct = /a/(l + fz)- For 
each redshift bin, we calculate the average logio Lx, and 
renormalise /ct to logioZ/x=43.5 in each bin. We present 
the data for these three redshift bins in Table [3] We plot 
the normalised Compton thick fraction as a function of red- 
shift in Fig. |14[ where we have merged the 1-2.5 and 2.5- 
4 redshift bins into one in order to improve the statistics, 
using a simple average. We also plot the results from Bur- 



lon et al.| ( |2011[ ) for the Swift/BAT sample of low redshift 
sources. The mean 2-8 keV luminosity for this sample is 
logiol/x/ergs _1 =43.6, once we have converted the 15-55 
keV luminosity using a F = 1.9 power-law, and is thus also 
comparable to the luminosity considered in our analysis. Our 
data point toward an increasing intrinsic CT fraction from 
«20% at z< 0.1 to ^40% at z> 1 for L x = 10 43 ' 5 ergs -1 
AGN. 

We also present number count data as a function of lim- 
iti ng flux (log N - l og S) for the full band selected sample 



of |Luo et al.| ( |2008| for AGN and CTAGN with L x > 10 42 
ergs -1 . These have been calculated as a function of 0.5- 
8 keV count rate, and then converted into a flux assum- 
ing a power-law model with T — 1.4 for convenience. The 
CTAGN numbers have been corrected for spectral identifi- 
cation incompleteness and contamination as described pre- 
viously, but the luminosity ranges have not been restricted 
as above, only to Lx> 10 42 ergs" 1 . Using these data, we 
have calculated the CT fraction as a function of limiting 
flux, which is very useful for comparing to XRB synthesis 
models. These data are presented in Fig. |15| along with the 
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Figure 15. Top panel: log N - log S data for AGN (solid line, 
where the dotted lines show the Poisson upper and lower uncer- 
tainties) and CTAGN (data points) with L x > 10 42 ergs -1 in 
the full band selected sample of Luo, et al (2008). Bottom panel: 
The Compton thick fraction as a function of limiting flux in the 
0.5-8 keV band. Fluxes have been converted from the 0.5-8 keV 
count rate using a V = 1.4 power-law. The dashed line shows the 
model prediction from Gilli, et al (2007). 



model prediction of Gilli et al. (20071. This fluxes from this 



prediction have been converted from a 2-10 keV flux to a 0.5- 
8 keV flux using a power-law model with T — 1.4. While this 
comparison is not straightforward, due to differing spectral 
models used, the observational and model data agree very 
well. However, the redshift evolution of the absorbed AGN 
fraction that we find here is not included in the XRB mod- 



els of Gilli et al. (20071, so the agreement we find with their 



models is likely to be coincidental. 



4-3.3 Alternative Compton thick selection methods 

A possible alternative method of identifying Compton thick 
sources, which as we have shown, can be very challenging 
from X-ray spectral analysis at low fluxes, is to use data 
from other wavelengths. Most studies have concentrated on 
mid-infrared (MIR) wavelengths, where the primary radia- 
tion absorbed by the circum-nuclear material is reemitted. 
Here, we briefly explore the multiwavelength properties of 
our sample starting with the MIR. We obtained 24 micron 
photometry in the CDFS from the FIREWORKS catalogue 



of Wuyts et al. (20081. We searched for counterparts to the 



X-ray sources by cross correlating the positions of the X-ray 
sources with the positions of the 24 micron sources, tak- 
ing the closest 24 micron source, rejecting any counterpart 
with an offset of > 5 arcsecs. Fig. [16] plots the X-ray flux 
against the 24 micron flux of X-ray sources with a 24 micron 
counterpart, highlighting broad line AGN, and the Comp- 
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Table 3. Details of the Compton thick fraction for AGN in the redshift range given by column (1) with an average 2-10 keV logioLx 
for all AGN given in column (2); Column (3) gives the total number of sources for each bin; Column (4) gives the number of A r n> 10 24 
cm -2 AGN in each bin identified from spectral fitting; Column (5) gives the number of A r pj> 10 24 cm -2 AGN in each bin after correcting 
for misidentification and contamination; Column (6) and (7) give the number density of AGN and CTAGN given the survey solid angle; 
Column (8) gives the Compton thick fraction in the selected luminosity range; Column (9) gives Compton thick fraction normalised at 
logioLx=43.5. The uncertainties presented are calculated from Poisson uncertainties. The uncertainties on /ct have been propagated 
from the Poisson errors on N(Nn> 10 24 cm -2 ) and N(Nu< 10 24 cm -2 ) 
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Figure 16. The 0.5-8 keV X-ray flux against 24 micron flux for 
all X-ray sources with a 24 micron counterpart. Blue squares are 
optical broad line AGN, and red squares are the Compton thick 
AGN we find here. The largest red squares are Compton thick 
quasars with an intrinsic Lx> 10 44 ergs -1 . 



Figure 17. The selection method of Fiore, et al (2008) for select- 
ing Compton thick candidates given red R-K colours and a MIR 
excess (selection criteria shown as a dotted box). Blue squares are 
optical broad line AGN, and red squares are the Compton thick 
AGN we find here. The largest red squares are Compton thick 
quasars with an intrinsic Lx> 10 44 ergs -1 . 



ton thick sources we find here. However, no clear separation 
of these sources can been seen. 

Secondly we explore the selection method of Fiore ct al. 



(2008J), whereby red sources in R-K (R-K > 4.5) with a MIR 
excess with respect to the optical (J 7 24 Hm /FR > 1000) are 
identified as Compton thick candidates. We take R band 



WFI photometry from Arnouts et al. (2001) and Ks band 
SOFI photometry from |Qlsen et al.| ( 2006] ~ Fig. [r7| shows 
these selection criteria and where the sources we study 
here lie. This method has been applied to undetected X-ray 
sources, and while only one source would be selected as be- 
ing Compton thick using this method, it is indeed Compton 
thick as shown in our analysis. 

Finally, the use of hardness ratios has often been used 



to infer the level of obscuration occurring when there are 
not enough counts for spectral fitting. These are computed 
using the number of counts detected in the hard band (H, 
typically 2-10 keV) and in the soft band (S, typically 0.5-2 
keV), the hardness ratio then being H-S/H+S. This is of 
course very crude and does not account for spectral com- 
plexity, such as soft excess emission. Brightman & Nandra 
(2012) have shown that a pair of rest frame hardness ratios, 
HR1 (2-4/1-2 keV) and HR2 (4-16/2-4), do very well at 
picking out very obscured sources (A?h> 10 23 cm -2 ) which 
are hard in HR2 due to absorption, and soft in HR1 due to 
the soft excess emission, which we have found to be ubiq- 
uitous here in this analysis. Figure [18] shows this selection 
method for the CDFS sources with 0.8 < z < 1.2, where 
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Noguchi et"aT| (|20To|> find that 25% of obscured AGN in 



Figure 18. The selection method of Brightman & Nandra (2012), 
which uses HR1 (F(2-4 keV)-F(l-2 keV)/F(2-4 keV)+F(l-2 keV)) 
and HR2 (F(4-16 keV)-F(2-4 keV)/F(4-16 keV)+F(2-4 keV)) to 
select heavily obscured sources, where the box has been shown to 
include all Compton thick AGN. Shown for CDFS sources with 
0.8 < z < 1.2. Blue squares are optical broad line AGN, and red 
squares are the Compton thick AGN we find here. The largest red 
squares are Compton thick quasars with an intrinsic Lx> 10 44 
ergs -1 . 



all the bands used in this scheme are observed in the Chan- 
dra band. The solid wedge shows where all Compton thick 
sources are expected to lie. The separation of Compton thick 
sources and BLAGN here is much cleaner than with the 
multi- wavelength approaches. 



5 DISCUSSION 

5.1 Geometrically buried AGN 

We find a significant number (~ 20%) of heavily absorbed 
sources in the CDFS in which evidence suggests that the 
central engine is buried in a close to An steradians of ma- 
terial. These sources are best fit by a torus model with a 
0° opening angle. Such geometrically buried AGN were pre- 
dicted to be a large fraction of the local AGN population by 
Ueda et al. (20071 from the discovery of two such sources 
in a then, very small Swift/BAT sample. They reach this 
conclusion from the measurement of scattered nuclear light, 
which is believed to be scattered by gas which fills the cone of 
the torus, being a very small fraction of the intrinsic emis- 
sion (/ S catt < 0.5%). Due to the higher redshift nature of 
our sources, we cannot place such strict upper limits on the 
scattered fraction as the soft X-rays get redshifted out of the 
Chandra band. However we do find five sources in which we 
can place an upper limit of / sca tt to be less than 0.5%. 



Earlier work by Levenson et al. ( 2002 I suggested the 



existence of AGN with tori that have small opening angles, 
using the equivalent width of the Fe Ka as a diagnostic, but 
degeneracies in the parameters meant they could not confirm 
this assertion. Their sources were found to be in composite 



starburst/AGN sources, supporting the proposal by Fabian 
et al. ( 1998 1 which suggests that starbursts in the inner 100 



pc of AGN are responsible for the nuclear obscuration in 



the XMM-Newton Serendipitous Source Catalogue are geo- 
metrically buried. Further evidence for absorbed AGN with 
very low / sca tt values comes from Eguchi et al. ( |2009[ ) with 
Suzaku observations of six local AGN and IComastri et aD 



I 2010) with Suzaku observations of five local AGN. Comastri 



et al. (20101 discuss the implications of these low / sca tt val- 



ues with respect to the XRB, concluding that the fits to the 
XRB spectrum and source counts are essentially unaffected. 
Other works that have found evidence for buried AGN in- 



clude Rowan- Robinson et al. ( 2009 1 who suggest candidates 



for AGN with close to 100% covering fractions in the Lock- 
man Hole from sources with high torus dust luminosities 
with respect to their X-ray luminosities. 

In this work, we attribute all soft excess emission in 
obscured sources to scattered nuclear X-rays. However, soft 
excess emission in obscured AGN can also be attributed to 
emission by circumnuclear gas photoionised by the AGN 



(Guainazzi & Bianchi 20071, or to star-formation related 



processes in the host galaxy, such as X-ray binaries. Due 
to the spectral quality of our data we are unlikely to be 
able to discriminate between the different emission compo- 
nents. Nevertheless our conclusion from / sca tt comes from its 
variation with ]Vh which is unlikely to be affected by these 
contaminants as photoionised emission is also likely to de- 
pend on the covering factor of the torus. While there have 
been suggestions that star formation and AGN obscuration 
are linked, the expected result from this is that /scatt would 
increase with ]Vh, rather than decrease. 

Our data show that the more obscured along the line 
of sight an AGN is, the more likely it is to be buried in 
a geometrically thick torus, shown by small scattered frac- 
tion in these sources, and the preference to be fit by torus 
models with small opening angles. A similar conclusion was 
reached by |Ricci et al.| ( ]2011[ ) , who found from the study 
of local Seyfcrt galaxies using hard X-ray data from INTE- 
GRAL, that Seyfert 2s with the strongest absorption also 
show the strongest reflection with respect to the less ab- 
sorbed Seyferts, which they attribute to a larger covering 
fraction of the torus in the most obscured sources. Our data 
does not however show that geometrically buried AGN are 
more prevalent at lower luminosities, or that the opening 
angle of the torus is dependant on luminosity. This would 
be expected when considering the 'receding torus' model 
( Lawrence||1991 Simpson]|2005 \ , which predicts that as lu- 
minosity increases, the covering factor of the torus decreases 
as luminosity increases due to the radius of dust sublima- 



tion increasing and has much observational evidence (Ueda 



et al. | 2003 . Hasinger 2008] |Burlon et al.|2011| |Brightman~ 



Nandra||2011b[ ). However, with the limited number of torus 
fits we have at hand, it may not be possible to break all 
degeneracies between Nh, £x and redshift. 



As noted by Elitzur ( 2012 1, the variation in torus cover- 



ing factors observed in AGN has important implications for 
the AGN unification scheme. This scheme employs only the 
orientation angle of the observer to explain the differences 



seen in AGN types ( Antonucci 1993 Urry & Padovani 1995 1, 



low luminosity AGN. Studies by Noguchi et al. (20091 and 



however it is now clear that the covering factor of the torus 
is also an important parameter. Open issues with AGN ob- 
scuration and unification are also reviewed in Bianchi et al. 
(2012), including the discussion of the now apparent need for 
a model to explain the complex phenomenology observed in 
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AGN. We have shown here that X-ray spectral analysis can 
give valuable information about the circumnuclear material 
in AGN, such as the covering factor, when the appropri- 
ate spectral models, such as those utilised in this paper, are 
exploited, and should be used in future to aid in a better 
understanding of AGN unification. 

5.2 Compton thick AGN 

Of the 20 secure Compton thick sources in the CDFS we 
find in this investigation, source numbers 137, 178, 180, 
186 and 282 were identified as Compton thick candidates 
by |Tozzi et ai] ((2006} and 96, 137 and 186 were identi- 



fied as Compton thick candidates by [Georgantopoulos et al.| 
(20071, both with lMs Chandra data. Furthermore, source 



number 309 was presented as a Compton thick source by 



Feruglio et al. 



(20111 using the 4Ms data. Meanwhile, Co 



mastri et al. (20111 presented two Compton thick sources 



from their analysis of ultra-deep XMM-Newton data in the 
CDFS, which are source numbers 176 and 265 in the |Luo| 



et al. (20081 catalogue. From our spectral fits to the Chan- 



dra spectra we find that for 176 Nn=6.57tiM x 10 23 cm" 



and 265 iV H =1.04 



+0.37 



< 10 J4 cm" 2 . Our result for 265 is 
to' x 10 24 cm" 2 ) to the fit to 
the XMM-Newton spectrum of the same source with the 



in very good agreement (1.0 



torus model of Murphy & Yaqoob ( 2009 1 by Comastri et al 



(20111. Of the remaining 21 l-a sources, 262, 265, 313 and 



337 were identified as Compton thick candidates by |Tozzi| 
et al. (20061 and 262, 313 and 337 were identified as Comp- 



ton thick candidates by Georgantopoulos et al. (12007). Thus 
we are presenting 29 sources here as Compton thick which 
have to our knowledge, not been identified through X-ray 
spectroscopy previously. 

We have accounted for the fraction of Compton thick 
sources which have not been identified through our spectral 
fitting, along with contaminating sources, with data from 
simulations. Furthermore, we take into account the selec- 
tion bias against Compton thick sources due to extreme flux 
suppression in the X-ray band, which we do by restricting 
the intrinsic luminosity ranges considered to those which 
we know the survey is sensitive to Compton thick AGN. 
We find that the Compton thick fraction in the CDFS is 
a moderately increasing function of redshift. The fact that 
the Compton thick fraction does depend on redshift is con- 



sistent with the suggestion from Treister et al.| (2009) that 
Compton thick sources have a steeper evolution from high 
redshift than less obscured sources. However, these authors 
report a much stronger increase in the Compton thick frac- 
tion of luminous AGN at z > 2 based on X-ray stacking 
analysis of infrared excess sources. Our conclusion is also 
consistent with results that find the obscured fraction of 
AGN (A r H> 10 22 cm -2 ) increases with redshift (|La Franca 



et al.||2005| |Ballantyne et al.||2006| |Treister fc Urry||2006J 
Hasinger|2008 1 . It is als o quantitatively consistent with the 



(z + l) u - az relation that |Hasinger| (2008) finds for f 2 vs red- 
shift, given the relation between f 2 and /ct that we have 
assumed. The Compton thick fraction in the local universe 
(z < 0.1) has been found to be «20% ( jBrightman fc Nandra| 
2011b Burlon et al.pOTl i from both mid-infrared selected 
samples and hard X-ray Swift/ BAT samples and our results 
show that this increases to ~40% at z=l-4. Our analyses in- 
clude sources with Ah> 10 25 cm -2 , which Chandra detects 



<8 keV through the scattered soft X-ray emission or re- 
flected emission. Swift/BAT is less sensitive to these sources 
due to extreme flux suppression in the >10 keV band. How- 
ever the results from the MIR selected sample of Brightman 
& Nandra (2011b I, which should be less biased against these 



extremely obscured sources, corroborate the hard X-ray de- 
tected results. 

The normalised intrinsic Compton thick fractions that 
we have calculated are dependant on the /ct vs. Lx relation 
that we have derived from the f 2 vs. Lx relation presented 
in Hasinger (20081. The assumption that we have made in 



deriving this relation is that the number of sources with 
10 24 <iV H < 10 26 cm" 2 is equal to that with 10 21 ' 5 <A H < 
10 24 cm -2 . If we were to change this and assume that the 
number of sources with 10 24 <Ah< 10 26 cm -2 is only half 
that of 10 21 ' 5 <A r H< 10 24 cm" 2 , then our normalised Comp- 
ton thick fractions would change by less that 4%, which is 
much less than the statistical error. 

Several authors have presented results on work to con- 
strain the space density of Compton thick AGN at high red- 
shifts, particularly at z«2, many of which include stacking of 
X-ray undetected sources. Of these works, our results com- 



pare well with those of Fiore et al. ( 2008 1, who find that the 



space density of infrared selected sources is in agreement 
with the XRB predictions of GCH07. Analyses by |Daddi| 
et al. ( 2007[ ) have on the other hand over predicted the space 
density with respect to the XRB models. |Alexander et al7| 
(20111 have placed X-ray spectral constraints on X-ray de- 



tected z~2 BzK selected galaxies in the CDFS and find that 
the Compton thick space density at these redshifts is lower 
than that predicted by the XRB. 

We also briefly investigate alternative selection meth- 
ods of Compton thick AGN, including multiwavelength se- 
lection. We find that Compton thick sources cannot be 
cleanly separated from the non Compton thick sources on a 
Fx — i*24|xm diagram as they occupy the same regions of the 



diagram as BLAGN. Georgantopoulos et al. (20111 explore 
the use of the Lx — ^6nm ratio as an effective way of select- 
ing Compton thick sources, and came to a similar conclusion, 
that this method cannot provide complete Compton thick 
AGN samples. While the selection method of |Fiore et al.| 
(2008) is aimed at X-ray undetected sources, it is clearly 
not suitable for sources detected in X-ray as only one source 



in the 2Ms catalogue of Luo et al. ( 2008 1 is selected as be- 



ing Compton thick by their selection criteria. A scheme by 



Brightman & Nandra (20121 based on hardness ratios of- 



fers a promising alternative to spectral fitting for selecting 
heavily obscured sources. 



6 CONCLUSIONS 

We have applied new X-ray spectral models which account 
for Compton scattering and the geometry of the circumnu- 



clear material from Brightman & Nandra (2011a I to spectra 



of high redshift AGN in the Chandra Deep Field-South. Our 
conclusions from this investigation are as follows: 

• We find evidence that the covering factor of the obscur- 
ing material and the line of sight Ah are linked, hence the 
most heavily obscured sources are also the most geometri- 
cally buried. We find this from the decline of the scattered 
fraction, / sca tt with Ah, and from the prevalence of the most 
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heavily obscured sources to be best fit by the torus model 
with the smallest opening angle. 

• We find that a significant number of heavily obscured 
(iV H > 10 23 cm~ 2 ) sources (~ 20%) in the CDFS are buried 
in geometrically thick material, having been fit best with 
the torus model with a 0° opening angle. 

• The average scattered fraction for sources in the CDFS 
is 1.7 ±0.4%. 

• We find a total of 41 Compton thick AGN in the CDFS 
from our analysis, 29 of which have not been identified as 
Compton thick by previous work on X-ray spectral analysis. 

• Binning our sample by redshift, considering only 
sources with an intrinsic luminosity greater than that to 
which the CDFS is sensitive, we find that the intrinsic 
Compton thick fraction, normalised to Lx= 10 43 ' 5 ergs -1 
increases from ss 20% in the local universe to « 40% at 
redshift 1-4. 
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